Universal vortex formation in rotating traps with bosons and fermions 
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When a system consisting of many interacting particles is set rotating, it may form vortices. 
This is familiar to us from every-day life: you can observe vortices while stirring your coffee or 
watching a hurricane. In the world of quantum mechanics, famous examples of vortices are super- 
conducting films* and rotating bosonic 4 He or fermionic 3 He liquids-^. Vortices are also observed 
in rotating Bose-Einstein condensates in atomic traps-*-*- and are predicted to exisS- for paired 
fermionic atoms- 4 -. Here we show that the rotation of trapped particles with a repulsive interaction 
leads to a similar vortex formation, regardless of whether the particles are bosons or (unpaired) 
fermions. The exact, quantum mechanical many-particle wave function provides evidence that in 
fact, the mechanism of this vortex formation is the same for boson and fermion systems. 
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Let us now consider a number of identical particles with 
repulsive interparticle interactions confined in a harmonic 
trap under rotation. These particles could be electrons 
in a quantum dolpifi, positive or negative ions, or neu- 
tral atoms in boson or fermion condensates 11 . Though 
simple to describe, this quantum mechanical many-body 
problem is extremely complex and in general not solv- 
able exactly. Consequently, in rotating systems the for- 
mation of vortices and their mutual interaction is usu- 
ally described using a mean field approximation. In 
superconductors this is the Ginzburg-Landau method 1 . 
For Bose-Einstein condensates, one often applies the 
Gross-Pitaevskii equationii*i£. In this way, Butts and 
Rokhsar— found successive transitions between stable 
patterns of singly-quantised vortices, as the angular mo- 
mentum was increased. A single vortex appears when 
the angular momentum L is equal to the number of par- 
ticles N, two vortices appear at L ~ 1.75 AT and three 
vortices at L ~ 2.1N (see Refsi 13 i 14 i 15 i 16 ). For quantum 
dots in strong magnetic fields, the occurrence of vortices 
was very recently discussed by Saarikoski et almL. 

Based on the rigorous solution of the many-particle 
Hamiltonian, we show that striking similarities between 
the boson and fermion systems exist: the vortex forma- 
tion is indeed universal for both kinds of particles, and 
the many-particle configurations generating these vor- 
tices are the same. For a small number of particles, 
the many-body Hamilton operator can be diagonalised 
numerically. We use a single particle basis of Gaussian 
functions to span the Hilbert space. These Gaussians are 
eigenstates of the trap for radial quantum number n = 
and different single-particle angular momenta. These 
states dominate for large total angular momenta L. We 
only consider one spin state, i.e. bosons with zero spin 
or spin-polarised fermions. Numerical feasibility limits 
calculations to small particle numbers TV. However, the 
advantages of our approach as compared to mean field 
methods are that our solutions (i) are exact (up to nu- 
merical accuracy), (ii) maintain the circular symmetry 
and thus have a good angular momentum, and (iii) allow 
the direct, quantitative comparison between boson and 
fermion states and thus serve to uncover the origin of the 
vortices in small systems. 



The many-particle energies for rotating clouds of bosons 
or fermions are compared in Figs. 1 (bosons) and 2 
(fermions). The low- lying states are shown as a function 
of total angular momentum L. Following the tradition 
in nuclear physics, the line connecting the lowest states 
at fixed L is called^ the "yrast" line. (The word "yrast" 
originates from Swedish language and means "the most 
dizzy".) It is marked with a red line in Figs. 1 and 2. 
For bosons the two spectra appear almost identical, al- 
though one of them is calculated with a short-range con- 
tact interaction and the other with a long-range Coulomb 
interaction^. 

The comparison between the lowest energy states of the 
spectra for bosons and Coulomb-interacting fermions, as 
displayed in Figs. 1 and 2, reveals striking similarities: 
The yrast line has the same kinks and vortices can be 
found in both systems appearing at similar angular mo- 
menta, as we explain below. 

When studying the appearance of vortices in the boson 
or fermion densities, we should remember that in con- 
trast to mean field methods, for an exact calculation 
with good angular momentum the particle density has 
circular symmetry and thus does not display the inter- 
nal structure directly. To find the vortices we therefore 
would need to study pair-correlation functions. In the 
fermion case, however, this can be problematic due to 
the disturbance of the exchange hole. Alternatively, as 
done in the insets of Figs. 1 and 2, we break the cir- 
cular symmetry with a small perturbation of the form 
Vg(r, ip) — Vq cos(£tp) which has t minima around the cen- 
ter. The perturbation can only couple states which differ 
in angular momentum by i£. Since the lowest energy 
state for each angular momentum is also the most impor- 
tant in the perturbation expansion, we can estimate the 
effect of the perturbation by the mixture of three yrast 
states = 9 (L) + rf[^fo(L-£)+^f (L + £)], where L 
is the angular momentum and rj is the mixing parameter, 
which is of the order 0.1 or smaller. If, for instance, the 
state in question has a two- vortex structure, then this will 
appear in ^(L) as two distinct minima for t = 2 already 
at very small mixing ratio rj. The insets to Figs. 1 and 2 
are obtained in this way. In the boson case, for N=8 the 
perturbative densities show a single vortex for L = 13, 
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FIG. 1: Rotating bosons in a trap. The Figure shows the low- 
lying many-particle energies of N = 8 bosons interacting by a 
contact interaction (yellow, lower left) or eight charged bosons 
interacting by the long-range Coulomb repulsion (white, upper 
right), as a function of angular momentum L (in units of K). 
The read line, also called the "Yrast line" , connects the lowest 
states at fixed L. The insets show vortices in the perturbative 
densities (as explained in the text), occurring above ratios 
L/N = 1 for a single vortex and L/N = 1.75 for two vortices. 



corresponding to L/N=1.6, while a second vortex occurs 
at L = 14, that is L/N = 1.75. For fermions, the angu- 
lar momentum is shifted^ by N(N — l)/2 as explained 
below. In analogy to the boson case, for N = 8 fermions 
the single vortex still exists at L = 28 + 13, while two 
vortices are found at L = 28 + 14. 

This universality in the vortex formation can be under- 
stood by looking more in detail at the many-particle 
states of the rotating system. In the case of non- 
interacting bosons the many-particle ground state is 



(i) 



where the coordinates in two dimensions are expressed 



by complex numbers Zj 



It turns out that 



in the case of spin-polarised Fermions the corresponding 
"condensate" is the so-called maximum density droplet 
(MDD)21 
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which is a Slater determinant of the consecutive single 
particle states, filled from m — to m^N — 1 and thus 
has a non-zero angular momentum Lp = N(N — l)/2. 
The state "J/ p corresponds to the Laughlin wave function 
for the integer quantum Hall effect. 
For bosons with short range repulsive interaction, a single 
vortex can be formed by multiplying the boson ground 
state by a symmetric polynomiaii£*2£ P 1V = Yik( z k~ z o), 



FIG. 2: Rotating, spin-polarised charged fermions (for ex- 
ample, electrons) in a trap. The Figure shows the low-lying 
many-particle energies for N = 8. The similarity to the boson 
case is remarkable. If bosons have a vortex state at angular 
momentum Lb , the fermions can have the same vortex struc- 
ture at angular momentum Lb + N(N — l)/2. In the example 
shown here, for fermions a single vortex can be seen in the 
perturbative densities above (L — 28)/N = 1. The inset shows 
a single vortex at L — 41, and two vortices at L = 42, corre- 
sponding to (L - 28) /N = 1.75. 



where zq is the center of mass. If one multiplies the MDD 
with the same polynomial 23 , this gives a good approxima- 
tion for the exact single vortex state for charged fermions. 
By noticing that for a system with many particles, the 
center of mass can be put at the origin, z$ = 0, we can 
make an Ansatz for the state with n fixed vortices form- 
ing a ring around the origin: 
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f„y = l[(z n - ae'" 1 ) x ■ • • x ]J(z jn - ae la ")^ B ,F 
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= n>? - a n )*B, F , (3) 



where ^ b,f is either the boson condensate or the fermion 
MDD, and the vortex centers are localised on a ring of 
radius a (oik = k ■ — ). This state does not have a good 
angular momentum, but such a state can be projected 
out by collecting only terms corresponding to a specific 
power of the constant a, 
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where S means symmetrisation. Note that with n = 1, 
K = N and a — this state describes a single vortex fixed 
at the origin. Figure 3 shows schematically the single 
particle occupation of these "vortex generating states" 
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for bosons and fermions. In both cases the n vortices are 
generated by exciting K single particles by n units of an- 
gular momentum. The quantum states of the numerical 
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FIG. 3: Vortex generating single-particle configurations for 
bosons and fermions. (The particle number is chosen to be 
eight in this example). Exactly the same excitations of the 
"ground states" (condensate for bosons, and maximum den- 
sity droplet for fermions) cause the vortices in both cases, 
with the same increase in angular momentum. 



N — 1 particles and plotting the phase as a function of 
the last coordinate 17 . This is done in Fig. 4 for the vor- 
tex generating configuration for N — 8. The state 
has maximum amplitude, when the electrons are located 
so that one electron is in the center and seven electrons 
form a ring around it. To study the phase, we fix six of 
the electrons on the ring and one at a slightly off-center 
position in the middle. The resulting phase is shown in 
Fig. 4a. One can clearly see that each electron carries a 
vortex with it, as known from the theory of the integer 
quantum Hall effect^. When the wave function is mul- 
tiplied with the polynomial generating the vortices, Eq. 
(4) with n = 2, two additional vortices appear (Fig. 4b). 
When the fermion state is replaced with the corre- 
sponding boson state <3> b only the two additional vortices 
are seen (Fig. 4c). 
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FIG. 4: The phase of the many-particle wave function for 
eight particles. The phase is shown as a function of the coor- 
dinates of one particle with all other coordinates fixed. The 
colour scale is such that the jump from — n to n corresponds 
to a change from orange to pink, (a) shows the phase of the 
fermion state of the maximum density droplet ^f, Eq. (2), 
with vortices localised at the fixed electrons, and (b) shows 
the phase of the two-vortex state. The two additional vor- 
tices appear in ^2V, when the state in (a) is multiplied by 
the vortex generating polynomial, see Eq. (4), and (c) shows 
the same two vortices for the boson state. 



exact solutions show that the dominating configurations 
in cases where we see 1, 2, or 3 vortices (independent of 
the number of particles) are indeed those shown in Fig. 3, 
for both bosons and fermions. 

In the exact diagonalisation, by multiplying the exact 
boson wave function with the wave function of the 
MDD, Eq. (2), we can determine the overlap between the 
fermion and boson states. It turns out to be even larger 
than the weights of the most important configurations. 
For example, the overlap between the two-vortex states 
shown in Figs. 1 and 2 is 57% while the weight of the 
most important configuration (as shown in Fig. 3) in the 
boson case is only 15% and in the fermion case 47%. 
The vortices are born by the rotational motion and con- 
sequently carry angular momentum. In the single par- 
ticle picture the angular momentum is associated with 
the phase of the complex wave function: Going around 
the angular momentum axis the phase changes by 2tt. 
Similarly, the phase changes by 27r in going around a 
vortex center. In the many-particle picture the phase of 
the wave function depends on the coordinates of all the 
particles. In this case, the phase change around the vor- 
tex cores can be visualised by fixing the coordinates of 



Finally, we will return to the fermion spectrum shown 
in Fig. 2. The maximum amplitude of the MDD, i.e. 
the fermion "condensate" , corresponds to the equilibrium 
particle positions of a classical system with logarithmic 
repulsive interactions^. In small systems a rigid rotation 
of this localised state gives some of the high angular mo- 
mentum statesi2i22iS&. For example, in the case of eight 
particles there are two classical configurations: a single 
ring of eight particles, which we label by (0,8), and a 
ring of seven particles with one particle at the center, 
(1,7). The former allows rigid rotation at angular mo- 
menta 28 + 8 = 36, 28 + 16 = 44, etc., while the latter 
at 28 + 7 = 35, 28 + 14 = 42, etc. (as marked by arrows 
in Fig. 2). In the fermion systems, in some cases both 
localisation and vortex structure coincide. This is for ex- 
ample the case at L — 42 for eight fermions, which in 
the pair correlation function shows two vortices as well 
as very weakly localised particles arranged in a (1,7)- 
configuration. 

To summarise, we have shown with exact solutions of the 
many-particle systems that the vortex formation in rotat- 
ing traps of bosons and fermions have universal features: 
(i) They appear at certain angular momenta determined 
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by the number of particles and number of vortices, (ii) formation does not depend on the shape of the repulsive 
the many-particle excitations generating the vortices are interaction between the particles, 
the same for bosons and fermions, and (iii) the vortex 
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